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Abstract 

Following a long tradition of physicists who have noticed that the Ising model 
provides a general background to build realistic models of social interactions, we 
study a model of financial price dynamics resulting from the collective aggregate 
decisions of agents. This model incorporates imitation, the impact of external news 
and private information. It has the structure of a dynamical Ising model in which 
agents have two opinions (buy or sell) with coupling coefficients which evolve in time 
with a memory of how past news have explained realized market returns. We study 
two versions of the model, which differ on how the agents interpret the predictive 
power of news. We show that the stylized facts of financial markets are reproduced 
only when agents are over-confident and mis-attribute the success of news to pre- 
dict return to herding effects, thereby providing positive feedbacks leading to the 
model functioning close to the critical point. Our model exhibits a rich multifractal 
structure characterized by a continuous spectrum of exponents of the power law 
relaxation of endogenous bursts of volatility, in good agreement with previous an- 
alytical predictions obtained with the multifractal random walk model and with 
empirical facts. 
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1 Introduction 



Many works borrow concepts from the theory of the Ising models and of 
phase transitions to model social interactions and organization (see, e.g., 
Callen and Shapero, 1974; Montroll and Badger, 1974). In particular, Orlean 
(1984, 1986, 1989a,b, 1991, 1995) has captured the paradox of combining ra- 
tional and imitative behavior under the name "mimetic rationality," by de- 
veloping models of mimetic contagion of investors in the stock markets which 
are based on irreversible processes of opinion forming. As recalled in the Ap- 
pendix, the dynamical updating rules of the Ising model are obtained in a nat- 
ural way to describe the formation of decisions of boundedly rational agents 
(Roehner and Sornette, 2000). The Ising model is one of the simplest models 
describing the competition between the ordering force of imitation or conta- 
gion and the disordering impact of private information or idiosyncratic noise 
(McCoy and Wu, 1973). In the same class of minimal models of emergent so- 
cial behaviors is the model of herding based on percolation clusters proposed 
by Cont and Bouchaud (2000). 

Starting with a framework suggested by Blume (1993, 1995); Brock (1993); 
Durlauf (1991, 1993, 1994, 1996, 1997, 1999), Phan et al. (2004) summarize 
the formalism starting with different implementation of the agents' decision 
processes whose aggregation is inspired from statistical mechanics to account 
for social influence in individual decisions. Lux and Marchesi (1999, 2000); 
Brock and Hommes (1999); Kirman and Teyssiere (2002) have also developed 
related models in which agents' successful forecasts reinforce the forecasts. 
Such models have been found to generate swings in opinions, regime changes 
and long memory. An essential feature of these models is that agents are 
wrong for some of the time, but, whenever they are in the majority they are 
essentially right. Thus they are not systematically irrational (Kirman, 1997). 

Here, we study a model of interacting agents buying and selling a single fi- 
nancial asset based on an extension of the Ising model. The agents make their 
decision based upon the combination of three different information channels: 
mutual influences or imitation, external news and idiosyncratic judgements. 
Agents update their willingness to extract information from the other agents' 
behavior based on their assessment of how past news have explained market 
returns. Agents update their propensity to herding according to what degree 
the news have been successful in predicting returns. We distinguish between 
two possible updating rules: rational and irrational. In the rational version, 
agents decrease their propensity to imitate if news have been good predictors 
of returns in the recent past. In the irrational version, agents mis-attribute the 
recent predictive power of news to their collective action, leading to positive 
self-reinforcement of imitation. We show that the model can reproduce the 
major empirical stylized facts of financial stock markets only when the updat- 



3 



ing of the strength of imitation is irrational, providing a direct test and the 
evidence for the importance of misjudgement of agents biased toward herding. 

Section 2 specifies the model and compares it with previous related versions. 
Section 3 presents the results of exhaustive searches in the space of the major 
parameters of the two versions of the model. We describe in turn the distri- 
butions of returns at multiple time scales, the auto-correlations of the returns 
and of the volatility (absolute value of the returns) at different time scales, 
the multifractal properties of the structure functions of the absolute values of 
returns and their consequences in the characteristic relaxation of the volatility 
after bursts of endogenous versus exogenous origins. Section 4 concludes. 



2 Model of imitation versus news impact 

2. 1 Definition of the model 

We study the following model of N agents interacting within a network Af 
(taken here for simplicity as the set of nodes linked by nearest-neighbor bonds 
on the square lattice; this implies that an agent sitting at a node interacts 
directly only with her four neighbors). At each time step t, each agent % places 
a buy (si(t) = +1) or sell (si(t) = —1) order. Her decision Si(t) is determined 
by the following process 



where E[sj] (t) is the expectation formed by agent i on what will be the decision 
of agent j at the same time t. The left-hand-side and the right-hand-side 
are in principle evaluated simultaneously at the same time to capture the 
anticipation by a given agent % of the actions of the other agents which are 
going to determine the change of the market price from t — 1 to t. Indeed, we 
assume that the decisions Sj(i) are formed slightly before t, in the period from 
t — 1 to t, when the news G(t) has become available in the interval from t — 1 
to t, and are then converted into price at t by the market clearing process. In 
principle, the best strategy for agent % is indeed to base her action for the next 
investment period on her best guess of the present action of all other agents 
for the next investment period (see below). 

Expression (1) embodies three contributions to the decision making process 
of agent i: 

• Imitation through the term Kij{t)E[sj]{t). The kernel is the 



Si (t) = sign ]T K i3 {t)E[ Sj }{t) + <Ti(t)G(t) + €i(t) 
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relative propensity of the trader i to be contaminated by the sentiment 
of her friend j (coefficient of influence of j on i). In other words, K^it) 
quantifies the strength of the influence of agent j's expected decision on the 
decision of agent i, which evolves with time as we soon specify. Due to the 
heterogeneity of the traders, Kij ^ Kji, generally. The sum J2jeAf i s carried 
over all agents j who are in direct contact with agent %. 
Impact of external news through the term ai(t)G(t). G(t) quantifies 
the impact of the external news I(t) on the decision of agent i. We follow the 
specification of the artificial stock market model formulated by Goncalves 
(2003), and assume that I(t) follows a standard normal distribution and 



G(t) = 



1, if I(t) > , 

w (2) 

-1, if I(t) < . 



<7i(t) quantifies the relative impact at time t of the news' positive or negative 
outlook on the decision process of agent i. In other words, <7j is the relative 
sensitivity of agent's sentiment to the news. 
• Idiosyncratic judgement associated with private information. €i(t) 
embodies the idiosyncratic content of the decision of agent % accounting for 
the interpretation of her own private information. We take it as being nor- 
mally distributed around zero with a constant (homogeneous) standard de- 
viation equal to unity, without loss of generality, since the relative strength 
of the three terms is already captures by the units of -fQ/s and oVs. 

The market price is updated according to 

p(t) =p(t-l)exp[r(t)] , (3) 

so that p(t) and r(t) are known at the end of the interval from t — 1 to t. The 
return r(t) is determined according to 

r(t) = , (4) 

where iV is the number of traders in M and A measures the market depth 
or liquidity and is taken constant. In expression (4), the decisions Sj(t) are 
formed slightly before t, in the period from t — 1 to t and are then converted 
into price at t by the market clearing process. 

We account for the adaptive nature of agents and their learning abilities by 
updating the coefficient of influence of agent j on agent i according to the 
following rule: 

Kij{t) = bij + onKijit - 1) + 0r(t - l)G(t - 1) . (5) 
In this, we follow the large literature on the rationality of imitation and 
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of adaptive behavior when lacking sufficient information or when this infor- 
mation seems unreliable (see, e.g., Brock, 1993; Brock and LeBaron, 1996; 
Brock and Hommes, 1997, 1998, 1999; Hommes, 2001; LeBaron, 2000; Kir- 
man, 1991; Lux, 1995, 1998; Lux and Marchesi, 1999, 2000; Takayasu, et al., 
1992; Youssefmir, et al., 1998; Levy, et al., 1995; Sato and Takayasu, 1998; 
Levy, et al., 2000; Gaunersdorfer, 2000). The coefficients quantify the in- 
trinsic imitation influence of agent j on agent i in absence of other effects. For 
cti = j3 = 0, we recover a constant coefficient of influence, which derives from 
a simple argument of bounded rationality recalled in the Appendix with (23). 
The coefficient cti > (possibly different from one agent to another) quanti- 
fies the progressive loss of memory of past influences on the present. The last 
term with f3 ^ quantifies how agent % updates her propensity for imitation 
based on the role of the exogenous news G(t) in determining the sign and 
amplitude of the observed return in the preceding time period. This update 
depends upon whether the news G{t — 1) known between time t — 2 and t — 1 
has the same sign as the price variations from t — 2 to t — 1, i.e., the same sign 
as the return r(t — 1) defined by (3). In addition to the sign, the amplitude 
of the return is also taken into account in the updating rule of the coefficient 
Kij (t) quantifying the propensity to imitate: indeed, a small amplitude of the 
return has low psychological as well as financial consequences and should not 
count as much as a large amplitude of the return. The simplest specification 
is to take into account the impact of the amplitude of the return linearly in 
its size, hence the form (3r(t)G(t — 1) in (5). Stronger nonlinear dependence is 
probably more relevant (Pandey and Stauffer, 2000; Ide and Sornette, 2002; 
Sornette and Ide, 2002) but is not considered further here to keep the dis- 
cussion simple. Krawiecki, et al. (2002) have considered a simpler version in 
which each (t) is purely random and is constructed as the sum of two ran- 
dom noises, one which is common to all coupling coefficients and one which 
is specific to it. They are able to reproduce volatility clustering and a power 
law distribution of returns at a single fixed time scale. 

Gongalves (2003) considered this model (5) with a.i = 0, i.e., with no memory 
of past influence on present influence, which leads the time series of Kij(t) 
looking like a white-noise process. Our addition of the memory effect modifies 
this white-noise structure into a Ornstein-Uhlenbeck type noise, tending to 
a random- walk-like process for ctj = 1. With the new term ai ^ 0, K^t) 
keeps a memory of past successes that the news had on predicting the stock 
market moves over approximately 1/(1 — oci) time steps, thus characterizes 
the strength of the persistence of the links of agent i with other agents. 

The sign of the coefficient (3 is crucial. 

(1) For (3 < 0, agent i is less and less influenced by other agents, the better has 
been the success of the news in determining the direction and amplitude of 
the market return. This process is self-reinforcing since, as decreases, 
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the dominant term becomes <Ji(t)G(t), which further ensures that the 
news correctly predict the decision of agents and therefore the direction 
of the market move, thus decreasing further the coefficient of influence 
Kij. Reciprocally, agents tends to be more influenced by others when 
the news seems to incorrectly predict the direction of the market. The 
news being not reliable, the agents turn to other agents, believing that 
others may have useful information (see below and the Appendix for an 
elaboration of this argument). 

Taking a negative j3 corresponds to agents behaving according to stan- 
dard rational expectations with respect to the flow of external news. If the 
stock market is in agreement with the news most of the time, a rational 
agent would conclude that the impact of imitation, of herding, of trend 
following and of other endogenously generated positive feedbacks, is mi- 
nor and the news are the dominating factor. Indeed, standard economics 
views the stock market as a machine transforming news into prices and 
the market is efficient when all news have been correctly incorporated 
and are continuously incorporated into the market prices. In our frame- 
work, this situation arises when the strategy Sj(t) = sign [(?(£)] consisting 
in following the news is found at least as good as or better than (1). 
The agent would in this case decrease its propensity K to imitate, and 
continue to do so as long as the news are most of the time in agreement 
with the stock market moves. We can thus summarize this case (3 < as 
describing "boundedly rational" agents. 
(2) For /3 > 0, the more the news predict the direction of the market, the more 
the agents imitate other agents. In other words, there is a reinforcement of 
the influence between agents when the news and the stock market return 
match at the previous period. This is the "irrational" case where agents 
either mis-attribute the origin of the market moves to herding rather than 
to the impact of news, or misinterpret the exogenous character of news in 
terms of endogenous herding or infer that other agents will be following 
more eagerly as a group the direction given by the news. 

The regime /3 > may result from several mechanism. 
• Mutually-reinforcing optimism. When the market is rising (r(t — 
1) > 0) and the news are good (G(t — 1) = 1), the agents may seek each 
other in order to determine if the rise can be sustained: if the agent's 
neighbors are all bearish, then the agent interprets this as a sign that 
the rise cannot be sustained due to a lack of majority support; if the 
friends are bullish, the agent is encouraged to feel optimistic, in the 
sense that more good news may be on the way. It is thus not so much 
the price rise or fall that agents try to predict and reflect but rather 
the continuation of good news (or bad news) in the future and whether 
they expect or not the good news to continue. Since a rise in the market 
is the result of more people feeling optimistic that more good news are 
on the way, then the optimism is spreading like an epidemic and the 
market rises, reinforcing the influence coefficient (called propensity 
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to be influenced by the felling of others by Gongalves (2003)). 
• Overconridence. Another mechanism is the tendency for humans to 
exhibit overconfidence in their abilities, either individually or as a group. 
If they see that the stock market has moved in the same direction 
as the news indicated, they may conclude that the information pro- 
vided by the other agents has been valuable, since the stock market 
is supposed to follow the rule of the majority. Agents may thus be 
tempted to increase their imitation behavior as long as news and stock 
market continue to match. Or said differently, they attribute a larger 
value in the prediction of the news which they interpret as being in- 
fluenced by the majority opinion. Indeed, experiments committed by 
Darke and Freedman (1997) show that a lucky event can lead to over- 
confidence. Heath and Gonzalez (1995) have also compared decisions 
made alone to decisions made following interactions with others, and 
shown that, while interaction did not increase decision accuracy or 
meta-knowledge, subjects frequently showed stable or increasing confi- 
dence when they interacted with others, even with those who disagreed 
with them (see also Savoy and Beitel, 1997; Slater and Rouner, 1992; 
Roberts and Castore, 1972; Roby and Carterette, 1965). A possible in- 
terpretation is that the interaction serves the role of rationalizing the 
subjects' decisions rather than collecting valuable information. In the 
same spirit, Sieck and Yates (1997) have shown that exposition to others 
of the rationale behind decisions increase markedly subjects' confidence 
that their choices were appropriate. 

In summary, the model contains the following general ingredients: (i) the 
agents make decisions based on a combination of three ingredients: imitation, 
news and private information; they are boundedly rational; (ii) traders are 
heterogeneous (Ky and <7j); (iii) The propensity to imitate and herd is evolv- 
ing adaptively as an interpretation that the agents make of past successes of 
the news to predict the direction of the market. 



2.2 Specification of the updates of expectations of other agents' decisions 

The model is completely specified once the algorithm, used to construct how an 
agent estimates her expectation E[sj] (t) of other agents' decisions in expression 
(1), is given. Three possibilities can be considered. 

• E[sj](t) = Sj(t — 1): agent i expects a persistence of the other agents' de- 
cisions, similarly to a martingale condition; in absence of any information 
other than the past actions, the next predicted action is the last one. This 
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leads to transform (1) into the prescription 



Si(t) = sign 



E KiAfyAt 



ai(t)G{t) + €i(t) 



(6) 



Without the term o"j(t)G(t), this was the model adopted by Johansen, et 
al. (2000) to develop a theory of herding to explain financial bubbles as 
regimes of strong imitation between agents. Kaizokji et al. (2002) have also 
studied the dynamics of a stock market with heterogeneous agents in the 
framework of a recently proposed spin model for the emergence of bubbles 
and crashes. The Appendix recalls the derivation by Roehner and Sornette 
(2000) showing how this specification is a natural consequence of bounded 
rationality of agents. 

E[sj](t) = Sj(t): each agent % has access to the information of the action 
of other traders instantaneously and she can not do better than use it. 
Computationally, the solution E[sj](t) = Sj(t) is self-consistent since Si(t) 
depends on Sj(t) which itself depends on the former. Such self-consistent 
formulation can be solved by using an iterative algorithm as follows 



E[ Sj ](t) 



lim s\ 



[k] 



(*) 



(7) 



where 



s i (*) = sign 



E 



[k] 



(t) + <Ti(t)G(t) + €i(t) 



(8) 



with, for instance, the starting condition s\ (t) = Si(t — 1). 
Information cascade along a given path i = 1, 2, • • • ,N where the ordered 
sequence of indices % encodes the curvilinear abscissa along a path linking 
all agents, in the spirit of the information cascades (Welch, 1992; Bikhchan- 
dani, et al., 1992). In such a scheme, E[sj](£) = Sj(t — 1) if j > % and 
E[sj](i) = Sj(t) if j < i when considering agent i. In other words, we have 



E[s j ](t)=H(j-i)s j (t 



H(i-j) Sj (t) , 



(9) 



where H is the Heaviside function. In the literature on adaptive learning, 
this procedure is said to be backward-looking, because it does not model 
the other agents' future behavior given all past information, for instance, 
because agents are not able to form beliefs about the other agents' future 
choices. These rather naive precedures have been suggested by psychologists 
and animal behaviorists in the 1950s and have been tested in laboratory 
settings. 
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2.3 Restriction to model (Gongalves, 2003) with oti = corresponding to the 
absence of memory of the coefficients 's of imitation 



Our model is a straightforward generalization of the artificial stock market 
model formulated by Gongalves (2003), in which the coefficients of the 
influence of j on i are taken identical for all j's and are updated as follows: 



Kij(t) = Ki(t) =h + 0r(t - l)G(t - 1) . (10) 



The coefficients b^s capture the "natural" propensity of humans for imitation, 
which may vary from agent to agent. There are drawn at the beginning of the 
simulation from a uniform distribution between and some maximum positive 
value and remain fixed thereafter during the dynamics. The main difference 
between (10) and our specification (5) is the absence of the persistence or 
memory of Kij(t) on its past values in (10). 

Based on numerical simulations and synthesis, Gongalves (2003) argues that 
the model (1) using the implementation (9) corresponding to an informa- 
tion cascade along a specific path and with (10) reproduce all the important 
stylized facts of the stock market only for f3 > 0. This is interesting but 
actually not quite correct. Consider, for instance, the following parameters 
&max = 0.22 ~ 0.24, a max = 0.14 ~ 0.15 and CV = 0.8 ~ 0.9, which are 
recommended by Gongalves (2003) to reproduce the main stylized facts (fat 
tailed distributions of returns, clustering of volatility, bubbles, crashes). The 
problem is that the time series of returns generated with these parameters 
have unrealistic bimodal distributions of returns, as shown in Fig. 1. The ori- 
gin of this bimodal structure results from the fact that the model explores the 
ordered regime of the Ising model too often, corresponding roughly speaking 
to an average coupling constant [Kj) larger that the critical Ising value K c . In 
this gents form a crowd with a majority opinion which can flip with time 

between ±1 due to the feedback process of the news and return dynamics on 
the coupling coefficient. In other words, the two bumps of the distribution of 
returns are direct signatures of the existence of the two spontaneous ordered 
states of the Ising model when the coupling coefficient is above its critical 
value. The dynamics of decisions is thus characterized by more or less random 
flips of crowds of agents between the two opinions. The resulting stylized facts 
can not therefore be taken as genuine signatures of a realistic dynamics of 
agents' opinions. 
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3 Simulations and results of the model with memory in the dy- 
namics of influence coefficients 

The origin of the unrealistic bimodal distribution of returns shown in Fig. 1 
has to be found in the rather artificial property that the influence coefficients 
Kij(t)'s, which are updated according to (10), lose instantaneously from one 
time to the next the memory of their past values. As a consequence, the 
coefficients Kij(t) fluctuate with time approximately as a white noise with an 
amplitude controlled by that of the returns and this gives rise to abrupt shifts 
of the majority opinion as explained in the previous section. This property 
is unrealistic to model the real persistence of interactions between agents. 
Indeed, in the real world, people are connected through social networks that 
do not instantaneously reshuffle in response to external effects as expression 
(10) describe. Rather, social connections evolve slowly and exhibit significant 
persistence, as documented in numerous studies (see, for instance, Suitor, et 
al., 1997; Wellman, et al., 1997, and references therein). Networks of investors 
communicating their opinions and sentiment on the stock market are similarly 
persistent. 

We thus turn to model (1) implemented with (6) or along a specific information 
cascade (9), as described in section 2.2 We assume the existence of a memory 
in the dynamics of the influence coefficients, which are updated according to 
the simplified version of (5) given by 

Ki(t) = h + aKi(t - 1) + (3r(t - l)G(t - 1) . (11) 

where the 6,'s are uniformly distributed in (0, & max ) at the beginning of the 
simulation and remain frozen thereafter. The value 6j of agent % represents her 
idiosyncratic imitation tendency. Compared with (5), the memory parameter 
a is taken the same for all agents. We use a 50 x 50 lattice as the geometrical 
implementation of the social network, in which the agents are located at the 
nodes and each agent interacts only with her four nearest neighbors. The sen- 
sitivity <7j of agent % to the global news G(t) is uniformly distributed in the 
interval (0, <7 max ). The value <7j of agent i is again specific to her and quan- 
tifies her susceptibility to be influenced by external news. The idiosyncratic 
or private information term €i(t) of agent i is drawn at each time step from 
a normal distribution with zero mean and a standard deviation s £j j which is 
also different from one agent to another: s £j j is chosen at the beginning of the 
simulation (like 6j and <7j) to characterize agent i according to a value equal 
to the sum of a common constant CV and of a uniform random variable in 
the interval [0,0.1]. 

In our simulations, we fix A = 40 (which determines the scale of the returns to 
a value comparable to that of empirical observations) and a = 0.2. We have 
also investigated other values a = 0.4, 0.6, and 0.8 and obtain similar results. 
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We explore the properties of the model in the parameter space of 6 max , ffmax 
and CV. There is no loss of generality in fixing |/?| = 1 to explore the relative 
importance of the term (3r(t — l)G(t — 1), since the typical scale of the K^s is 
set by 6 max whose amplitude is varied in our numerical exploration. However, 
the sign of (3 is fundamental as explained in section 2.1. We thus consider in 
turn the two cases (3 — — 1 and (3 — +1 and explore in each case a large sample 
of triplets (b max , a max , CV). 

We ask whether the model can account for the most often reported stylized 
facts of financial markets. In other words, we would like to validate the model. 
For this, we consider the following metrics: (i) the distribution of returns 
at different time scales; (ii) the correlation function of returns and of the 
absolute value of the returns (taken as a proxy for the financial volatility); (iii) 
the scaling of the moments of increasing orders of the absolute values of the 
returns (testing multifractality); (iv) the existence of a hierarchy of exponents 
controlling the relaxation of the volatility after an endogenous shock (another 
hallmark of multifractality); (v) the existence of bubbles and crashes and their 
properties. Our strategy is to search for a robust set of model parameters 
for which all these stylized facts are reproduced not only qualitatively but 
also quantitatively. The main result of our analysis is that it is impossible to 
validate the model for (3 < while we find sets of parameters for (3 > which 
nicely fit the stylized facts of real financial markets. 



3.1 News predicts the next return — > decrease of imitation: (3 = — 1 

In this case, when the news G happen to correctly predict the return (rG > 
0), the agents reduce their mutual imitation. On the contrary, when rG < 
0, imitation among agents strengthens. Thus, (3 < corresponds to agents 
behaving according to standard rational expectations with respect to the flow 
of external news: a rational agent would concluded that, if the stock market 
is in agreement with the news then, the impact of imitation is minor and 
the news are the dominating factor, as standard economic textbook describe, 
since the strategy Si(t) = sign[G(i)] consisting in following the news is found 
as good as or better than the more elaborate dynamics of the agents' actions 
Si(t) incorporating the imitation process. For f3 < 0, agents decrease their 
propensity to imitation as long as the news are in agreement with the stock 
market moves. 

The following argument shows that the attractor of the dynamics is character- 
ized by negligible imitation and only the news and private information terms 
are important for the dynamics. Consider a population of traders at time t 
with their propensity Ki(t) to imitate on average above the critical Ising value 
K c such that imitation initially dominates the dynamics. The news G(t) be- 
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come known and the decisions of the agents given by (1) decide collectively 
the return r(t + 1) from time t to t + 1 through expression (4). Since the 
.£Q(i)'s are overall above K c , it is well-known from many past studies of the 
Ising model (see for instance McCoy and Wu (1973) and references therein) 
that the corresponding "ferromagnetic" phase of the system is characterized 
by a strong slaving to external fields such as G(t). Hence, the collective opin- 
ion J2i s i takes the sign of G(t) with a high probability. As a result, with a 
large probability, r(t + l)G(t) is positive, which entails a downgrading of Ki 
by an amount r(t+ 1) since beta = — 1 (in addition to the other terms which 
tend to reverse K^t) to the value 6^/(1 — a)). This behavior continues as long 
as the KiS are above or close to K c (since the K^s are heterogeneous, the 
effective critical values is modified compared with the homogeneous case and 
our argument remains valid when using this modified value). Alternatively, 
if the Ki's are on average smaller than K c , the collective decision J2i s i an d 
therefore the market return have little or no relationship with the external 
news. Hence, the term f3r(t + l)G(t) takes random signs from one time step 
to the next, leading to a effective random forcing added to the autoregressive 
equation K^t) = fej + aK^t — 1). The coefficients K^t) evolve to fluctuate 
around the asymptotic value — a). We thus expect Gaussian distributions 
of returns when — a) is smaller than K c and bimodal distributions when 
bi/(l — a) > K c reflecting the slaving of the global opinion to the sign of the 
news. 

In our simulations, we have scanned 6 max from 0.1 to 0.5 with spacing 0.1, cr max 
from 0.005 to 0.08 with spacing 0.005, and CV from 0.1 to 1.1 with spacing 0.2. 
This corresponds to a total of 480 different models. We use (1) implemented 
according to the information cascade of sentiment formation explained in the 
item associated with equation (9). This choice is made to minimize the com- 
putational cost. Tests using the other updates over limited time span suggest 
that the results we are interested in are not sensitive to the details of the 
updating rules. We run each model over 10 4 time steps and then analyze the 
time series of returns. 

The first metric we analyze is the distribution of returns at different time 
scales r, defined according to 

r T (t) = \n\p(t)/p(t - r)} . (12) 

where r is a multiple of the time step. We observe two classes of shapes. For 
large idiosyncratic noise (large CV) and not too large 6 max , the distribution of 
returns is Gaussian for all time scales r. For smaller CVs and larger 6 max , we 
observe multimodal return distributions, as illustrated by the typical example 
shown in Fig. 2 obtained for 6 max = 0.2, <j max = 0.045, and CV = 0.1. The 
number of peaks in the distribution of r T is r + 1. These multimodal distri- 
butions correspond to the regime where the news G(t) controls the collective 
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opinion of the traders which tend to coordinate their decisions as explained 
above when 6j/(l — a) is larger than K c . For given 5 max and CV, the bi- 
modal structure of the distribution of r\ becomes more and more significant 
as (T max increases. Alternatively, for fixed a max and CV (say, er max = 0.02 and 
CV = 0.1), a bimodal distribution is obtained for sufficiently small 6 max and 
we observe a crossover from bimodal (6 max = 0.1) to unimodal (6 max = 0.5) 
with a crossover with a plateau (6 max = 0.3). This corresponds to the regime 
in which a majority of agents react to the global news in the same manner and 
buy or sell simultaneously. Since the news are G(t) = ±1 with equal proba- 
bility, the decisions of the agents and therefore the returns most often jump 
between two values of equal amplitude and opposite signs. The multi-modal 
structure of the distributions of returns at time scale r then results from the 
properties of the convolution of the bimodal distribution of the returns r x . 

In the parameter space that we have explored and notwithstanding our best 
attempts, we have not been able to find a set of parameters leading to distri- 
butions of returns exhibiting a monomodal shape with fat tails for small time 
scales, evolving slowly towards Gaussian distributions at large time scales, as 
can be observed in empirical data (Ghashghaie, et al., 1996; Mantegna and Stan- 
ley, 2000). 

In addition, we observe that the correlation function of returns (C T (r,r)) and 
of volatilities (C T (|r|, |r|)) have similar amplitudes and decay with the same 
characteristic time scale as a function of time lag. This is very different from 
the observed correlations of financial markets, with very short memory for 
returns and long-memory for the volatility. 

We have also investigated the impact of the updating rule of the agents' sen- 
timents. If instead of the information cascade (third item of section 2.2), we 
use the parallel update (6), for the same range of parameters 6 max , cr max , and 
CV, we find that most of the returns are between two values proportional 
to ±1 and the distribution of returns is close to P{t\) = 5{r\ — 1). In other 
words, the bimodality of the distribution of returns is much more pronounced 
than shown in Fig. 2 that was obtained for the information cascade updating 
scheme. For large idiosyncratic noise (large CV), the bimodality disappears 
and is replaced by pdf's which are approximately Gaussian. With CV and 
&max (resp. <r max ) fixed, the volatility increases with increasing <7 max (resp. de- 
creasing 6 max ) and the distribution is bimodal when a max (resp. 6 max ) is large 
(resp. small) enough. The time evolution exhibits in addition long transients 
with returns fluctuating around zero before bifurcating to the bimodal state. 
The duration of the transient decreases with increasing <x max or 6 max . Scanning 
the parameters, we never obtained realizations with realistic distributions of 
returns at different time scales. We conclude that the updating scheme (11) 
with f3 — — 1 corresponding to boundedly rational agents cannot explain the 
stylized facts of empirical finance. 



14 



3.2 News predicts the next return — > increase of imitation: (3 = 1 

In our simulations, we fix a = 0.2, 6 max varies from 0.1 to 0.5 with spacing 
0.1, a max from 0.01 to 0.08 with spacing 0.01, and CV from 0.1 to 0.7 with 
spacing 0.2. This gives 160 models that we explore by generating time series 
of length equal to 10000 time steps. We use (1) implemented according to the 
information cascade of sentiment formation explained in the item associated 
with equation (9). 

When CV is very large, the distribution of return r T is Gaussian, simply be- 
cause the returns are dominated by the idiosyncratic noise modeling private 
information, which is chosen Gaussian. This regime is not interesting since it 
erases both the effect of news and of imitation. For smaller CVs, we also ob- 
serve bimodal distributions of the returns r 1 for certain ranges of parameters, 
as in the case f3 — — 1. This occurs for large a max and small 6 max . Similar bi- 
modal distributions of returns are obtained for a = 0, as described in section 
2.3. 

We have found several parameter combinations which lead to realistic stylized 
facts. For instance, for the following sets of (6 max , cr max , CV), (0.3, 0.03, 0.1), 
(0.4, 0.04, 0.1), (0.4, 0.05, 0.1), (0.5, 0.06, 0.1), (0.1, 0.01, 0.3), (0.1, 0.02, 0.3), 
(0.2, 0.02, 0.3), (0.2, 0.03, 0.3), (0.3, 0.04, 0.3), (0.5, 0.05, 0.3), (0.5, 0.07, 0.3), 
(0.3, 0.03, 0.5), and (0.5, 0.05, 0.5), the distribution of returns r T is a stretched 
exponential (or close to a power law) (Malevergne, et al., 2003) for small r, 
exponential for intermediate r, and Gaussian for large r, as in real financial 
data. 

In the following, we exemplify the obtained stylized facts with a = 0.2 and 
with the parameter combination (6 max = 0.3, (x max = 0.03, CV = 0.1) which is 
typical. 

3.2.1 Probability density functions of log-returns at different time scales 

Figure 3 shows a realization of the logarithm of the price over a time interval 
of 10 5 time steps. Figure 4 shows the corresponding time series of the log- 
returns defined by (3,4), where A = 40 is set to scale the returns to realistic 
values. Note the existence of clusters of volatility which are qualitatively sim- 
ilar to those observed in real financial data. The solid lines in Fig. 5 show the 
logarithm of the probability distribution densities (pdf) of the log-returns at 
different time scales defined by (12) as a function of the returns r T scaled by 
their standard deviation a T . The pdf curves have been translated vertically for 
clarity. In the semi-log representation of Figure 5, a straight line qualifies an 
exponential law. We observe stretched exponential laws at short time scales 
that cross over smoothly to a Gaussian law at the largest shown time scale. 
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This evolution of pdf's with time scales comply with the well-known stylized 
fact of financial markets (Ghashghaie, et al., 1996). The model thus obtains 
price series with the correct monomodal shape with fat tails, and the correct 
progressive transition to a Gaussian distribution at large time scales. We stress 
that this comparison improves on those involving only one time scale. 

It is interesting to compare the obtained pdf's at time scales r larger than the 
elementary time step 1 with those which would derive by r-fold convolution 
of the pdf at the unit time scale. In absence of time dependence, the two 
should be asymptotically identical. For instance, the pdf pdf(r 4 ) of returns 
r 4 over four time steps would be given by pdf(r 4 ) = pdf(ri) <g> pdf(ri) <S> 
pdf(ri)<E)pdf (r±). Upon convolution in the absence of dependence, the variance 
is additive which gives a prediction for the standard deviation of the pdf of 
returns at r time steps, a T = (Tiy/r, which can be used to normalize the pdfs 
in terms of the reduced variable r T /a T . The corresponding pdf's of r 4 , ri 6 , 
r 64 , and r 2 ^ obtained by convolution of the pdf of r\ are drawn in Fig. 5 
with dashed lines. We see clearly that the true pdf's for r = 4, 16, 64 have 
a much fatter tail compared with the theoretical pdf's based on absence of 
dependence. Such behavior is very similar to what is observed in real data 
(see for instance Malevergne and Sornette, 2005, Fig. 2.2). 

3.2.2 Autocorrelations of log-returns and volatility 

Figure 6 shows the temporal correlation of the log returns r\ as a function 
of the time lag £. One can observe a very short correlation time, of duration 
smaller than one time step. Figure 7 presents the temporal correlation of the 
absolute value of log returns r±, taken as a proxy for the volatility, both in 
linear-linear and in linear-log scales. It is apparent that the volatility exhibits 
a strong correlation with memory lasting approximately 100 time steps for this 
set of parameter. The right panel of Figure 7 represents the correlation of the 
volatility as a function of the logarithm of the time lag. This representation is 
suggested by the multifractal random walk (MRW) model which is constructed 
by definition with a correlation decaying linearly with the logarithm of the 
time lag, up to a so-called integral time scale T (Muzy, et al., 2000; Bacry, et 
al., 2001; Muzy, et al, 2001; Sornette, et al., 2003). The right panel of Figure 
7 shows that this dependence suggested by the MRW provides a reasonable 
approximation of the numerical data. The integral time scale is here estimated 
around 100 time steps. These observations are is good agreement with the 
stylized facts on the correlation of returns and of volatility of real financial 
markets. Indeed, one of the key stylized facts observed empirically is that there 
are only very short-range correlations in price changes and the time memory is 
less than one trading day and as small as minutes for the most liquid markets 
(Liu, et al., 1999; Gopikrishnan, et al., 1999; Mantegna and Stanley, 2000). In 
contrast, volatility exhibits a memory over up to the order of one year. 
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If we combine the information on the time scales from the pdf's of returns 
shown in Figure 5, the correlation of returns shown in Figure 6 and the cor- 
relation of the volatility shown in Figure 7, we can obtain a rough idea of 
the correspondence between the time step of the model and real trading time. 
From Figure 5, we see that the pdf of returns at time scale r = 4 and r = 16 
are similar to the empirical one at the daily scale for major stocks and indices. 
This suggests that one day corresponds to roughly 4 — 16 time steps of the 
model. This correspondence is compatible with Figure 6 for the absence of 
correlation at the daily time scale observed empirically. This correspondence 
gives with Figure 7 an integral time scale for the volatility correlation of about 
100/4 to 100/16 days, i.e., 6 — 25 days. This is about a factor of ten shorter 
than observed on real markets, if we believe the relevance of the MRW and 
its calibration of T for real markets. 

The memory of the autocorrelation of the volatility (as well as the correlation 
of the returns) is sensitive to the value taken by the parameter a introduced in 
our model, which embodies the dependence of the coefficient of imitation 
on its past values. Figure 8 shows that much longer ranges for the correlation 
of the absolute value of returns are found for larger values of a, however at 
the cost of introducing an unrealistic correlation of the returns. This figure 
suggests that a cannot be larger than 0.2 — 0.3 without producing unrealistic 
correlation in returns. 



3.2.3 Multifractal properties 

The MRW also predicts (and this is well-verified by empirical data) that the 
autocorrelation functions of |r T (t)| for different r should superimpose for time 
lags larger than their respective r (Muzy, et al., 2000; Bacry, et al., 2001). 
Figure 9 shows that this is approximately the case. 

Another important stylized facts is the multifractal structure of the abso- 
lute values of log-returns (Fisher, et al., 1997; Mandelbrot, 1997; Vande- 
walle and Ausloos, 1998; Brachet, et al., 2000; Bershadskii, 1999; Ivanova and Aus- 
loos, 1999; Schmitt, et al., 1999; Pasquini and Serva, 2000; Muzy, et al., 2001), 
which led to the proposition that the MRW might be a good model for finan- 
cial price time series (Muzy, et al., 2000; Bacry, et al., 2001; Muzy, et al., 2001; 
Sornette, et al., 2003). Figure 10 shows in log-log scale the structure function 

M q (r) = <|r r H (13) 

as a function of the time scale r. The power law dependence of the structure 
functions M q (r) as a function of r is found to be reasonable. The slopes of the 
lines in log-log plots give the exponents £ r defined by 

M q (r) ~ r*« . (14) 



17 



Multifractality is qualified by a nonlinear dependence of £ q as a function of the 
order q of the structure function (Mantegna and Stanley, 1995; Ghashghaie, et 
al., 1996)., as reported in Fig. 11. 

3.2.4 Endogenous versus exogenous shocks 

The dynamical process described by (5) together with (6) and (4) describes 
a flux of external news G(t) which are "digested" by the collective behav- 
ior of the population of traders to create a time series of returns presenting 
long-range memory in the volatility and multifractal properties, similarly to 
the multifractal random walk model (Muzy, et al., 2000; Bacry, et al., 2001; 
Muzy, et al., 2001; Sornette, et al., 2003). In addition, Sornette, et al. (2003) 
have discovered a new consequence of multifractality in the form of a contin- 
uous dependence of the exponent of the power law relaxation of the volatility 
after a spontaneous peak as a function of the amplitude of this peak (see also 
other applications in Sornette and Helmstetter, 2003; Sornette, et al., 2004). 
We proceed to test if such an effect is found in our model. 

Consider a realization giving a time series of returns r(t). Let us define the 
local volatility within a window [t + l,t + At] of size At by 

At 

^(*) = EK* + 0l 2 (is) 
i=i 

Using this definition, we construct a time series for the volatility o~\ t (t) in 
moving windows of size At. The (unconditional) average volatility E[a 2 ] is 
nothing but the average of the time series of o~\ t {t). Let us follow Sornette, et 
al. (2003) and consider local burst of volatility cr\ t {t) with amplitude scaled 
to the average volatility -E[o~ 2 ]: 

*l t ( t) = e 2 ^E[a 2 } . (16) 

The parameter s thus quantifies the relative amplitude of a local burst of 
volatility in units of the average volatility. Following Sornette, et al. (2003), 
for a given s, we identify all times t s whose volatility o~\ t {t s ) is close to e 2s E[a 2 }, 
that is, 

e 2{s ~ ds) E[a 2 ] < a 2 At (t s ) < e 2(s+ds) E[a 2 } , (17) 

where ds <^ 1. For a given relative log- amplitude s, we translate and superim- 
pose all time series starting at all the previously found times t s . Averaging over 
these time series of volatility obtains the average conditional relaxation func- 
tion of the volatility -E[cr(t|s) 2 ] following a local burst of volatility of amplitude 
(16). The multifractal random walk model predicts a power dependence 

E[a(t\s) 2 } ~ r aW , (18) 
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with 

" (S) = 3/2 + HT/At) ' (19) 

when At < t <^ Ate' s '/ A2 , where A 2 ~ 0.02. We keep the symbol a(s) for 
the exponent in (18) in line with the notation of Sornette, et al. (2003), but 
this should not be confused with the parameter a in (5) which controls the 
memory of the imitation coefficients K^. 

Fig. 12 shows the average normalized conditional volatility E[a{t\s) 2 ]/ E[a 2 ] 
as a function of the time t — t s to the local burst of volatility at time t s for 
different log-amplitudes s in double logarithmic coordinates. As expected from 
the predictions (18) and (19) of the MRW, the average volatility after a burst 
decays (respectively increases) when s > (respectively s < 0). Similar to 
real data analyzed by Sornette, et al. (2003), we observe approximate power 
laws. The power law exponents for different values of At are plotted in Fig. 13. 
The exponent a(s) depends linearly on s, with a slope increasing with At, in 
agreement with the prediction of the MRW and the finding of Sec. 3.2.3. 

We can actually obtain a direct estimation of the integral time scale T by 
studying how the slope 1/k of a(s) as a function of s depends upon At. 
Expression (19) predicts that k(At) = -±ln(At) + \ ln(T) + §. A linear re- 
gression of k as a function of ln(At) gives k(At) = — 0.481n(At) + 2.08. The 
first coefficient —0.48 is nicely close to the exact value —1/2 predicted by the 
multifractal theory, which provides an independent check on the validity of 
multifractality. Identifying \ ln(T) + § with 2.08 yields T = 14.3. 

Fig. 14 shows the average relaxation of the volatility after an exogenous shock, 
created by imposing a very large news impact G(t s ) at a single time t s and 
then letting the system evolve according to its normal dynamics thereafter. 
To gather sufficient statistics, we impose such large shocks with a periodicity 
of several hundred time steps, which is sufficiently long to allow the system 
to relax back to its normal fluctuating volatility. We have checked that the 
relaxation shown in Figure 14 is independent of the amplitude of the shock 
when sufficiently large. Note that, immediately after the news impact, the 
volatility first increases over a few time step before relaxing, which reflects 
the strong increase of the coupling coefficient Ky and the resulting stronger 
cooperativity of the agents. In order to mimic the previous analysis of the data 
by Sornette, et al. (2003), the origin of time for the relaxation of the volatility 
is taken at the peak time, rather than from the incipient exogenous news shock 
at t s . This shifts the origin of time by approximately 2 — 3 time steps, which 
is compatible with the correspondence that one trading day corresponds to 
4—16 time steps, as discussed in section 3.2.2. 

According to the theory relating the endogenous relaxation to the exogenous 
response function of the MRW developed by Sornette, et al. (2003), the re- 
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laxation shown in Fig. 14 should be characterized by an exponent close to 
1/2 over approximately the same range of times t — t s as found for the power 
law dependence of the relaxations shown in Fig. 12. We indeed observe a first 
decay regime which is compatible with a power law with an exponent close 
to 0.5 (but of course the range is too short to provide anything other than an 
indication). We also observe a cross-over to a faster decay, compatible with 
a faster decaying power with exponent close to 1.5. This behavior is actually 
expected if the system is not exactly critical but close to critical, as shown 
by Sornette and Helmstetter (2003) and Sornette, et al. (2004): the response 
function to an exogenous shock should in this case cross over from a depen- 
dence proportional to l/t 1 ^ 9 to l/t 1+e , with 9 = 1/2 for a multifractal system. 
Fig. 14 thus suggests that the multifractal properties of our system hold only 
up to a finite time scale T beyond which a cross-over to a non-critical behavior 
dominates. 



4 Discussion 

In this paper, we have extended the artificial stock market model introduced 
by Gongalves (2003) to include a memory in the dynamics of the influence 
coefficient on its past realization. This additional memory turns out to be a key 
ingredient to reproduce the major stylized facts of financial stock markets. In 
the previous specification a = of Goncalves (2003), the influence coefficients 
Ki adjust instantaneously to previous news and returns realization. With a 
non-zero a as proposed here, the influence coefficients exhibit an inertia. We 
believe that this is a crucial property of the interactions between social agents 
who only relatively slowly update their tendency to imitate their colleagues or 
friends. That this parameter provides, together with the competition between 
imitation and news, the main stylized facts of financial stock markets is an 
encouraging sign that we have correctly captured some of the most important 
ingredients at the origin of the organization of financial stock markets. 

These ingredients comprise imitation between agents, their influence by ex- 
ternal news and the impact of their private information. The imitation plus 
the idiosyncratic part of the decision process give together the dynamics of 
the Ising model. The news act then as an time-dependent external field. The 
addition of an evolution in the influence coefficients (which can also be called 
"coupling coefficients" ) make the strength of the imitations between agents a 
function of the past realization of the news and returns. 

The empirical stylized facts of financial stock markets have been found only 
for (3 > and a neither too small nor too large. The condition (3 > means 
that agents increase the propensity to imitate if the external news have been 
predictive of the returns in the past. This behavior corresponds to agents who 
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misinterpret, or misattribute the source of the prediction of returns. Alter- 
natively, this behavior corresponds to over-confident agents. Technically, the 
stylized facts in this regime result from the fact that the model operates around 
the critical point of the corresponding Ising model, with coupling coefficients 
which are time-dependent and endowed with a memory of past realizations. 
The critical point of the Ising model is associated with a critical value K c 
for the average coupling coefficient. Close to this value, agents organize spon- 
taneously within clusters of similar opinions, which become very susceptible 
to small external influences, such as a change of news. This may explain the 
occurrence of crashes, as argued previously (see, Sornette, 2003, Ch. 5). The 
present model exhibits bubbles and crashes as shown in figure 15, whose de- 
tailed study will be reported elsewhere. Intuitively, the critical slowing down 
well-known to characterize the proximity to the critical Ising point can explain 
the long-term memory of the volatility while the almost absence of correlation 
of the returns themselves is ensured by the impact of the news and the random 
idiosyncratic decisions. 

The fact that the most basic stylized facts (monomodal shape with fat tails, 
short-time return correlations and long-memory of the volatility) cannot be 
obtained for (3 < suggests the importance of the imitation behavior captured 
by (10) with /3 > 0. In this model, the creation of anomalous volatility, of 
its persistence, and of mult ifr act ality result from the tendency of agents to 
misinterpret the combined information of the news and of the stock market as 
resulting from the influence of the agents and their imitation. Or seen from a 
different view point, conditioned on their role of reflecting the stock market, 
the news serve as the substrate for fostering social interactions and reinforcing 
herding. By the mechanism of intermittent reinforcing of social interaction in 
(10), the coupling coefficient will vary and sometimes increase close to or 
cross a critical value at which critical fluctuations occur and beyond which 
global cooperativity dominates. 

As a bonus, we have discovered that this simple model exhibits a rich multifrac- 
tal structure, diagnosed not only by the standard convexity of the exponents 
of the structure functions but also by distinct power law response functions to 
endogenous compared with exogenous volatility shocks (Sornette, et al., 2003; 
Sornette and Helmstetter, 2003; Sornette, et al., 2004). To our knowledge, 
this is the first nonlinear model in which such clear distinction is documented 
quantitatively, based on a bottom-up self-organization. In contrast, the multi- 
fractal random walk which has provided the theoretical predictions used here 
is a descriptive phenomenological model. 

Let us end by noting the connection with the model of Wyart and Bouchaud 
(2003), which can be embedded in our model by putting all to (no imita- 
tion) and adding a dependence of Oi on past correlations between the realized 
returns r(t) and the available information G(t), similarly to the dynamics 
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of Kij in our model. Wyart and Bouchaud (2003) have studied the limit of 
self-fulfilling conventions created by the belief of agents on the existence of 
correlations between information and returns. By this belief, traders try to 
estimate this correlation from past time series and act on it, thus creating it. 
Our model emphasizes the other class of conventions based on imitation and 
moods. For the future, it would be interesting to combine both mechanisms 
as they are arguably present together in real markets, in order to clarify their 
relative importance and interplay. 



Appendix: Justification that the imitation term is rational in the 
absence of reliable information 

The form (6) of the formulation of the decision of agent % (without the term 
<Ji(t)G(t)) derives naturally from an argument of bounded-rationality, as fol- 
lows. 

Let us denote N(i) the number of traders directly connected to % on the graph 
of acquaintance (iVj = 4 for the 2D-square topology). The traders buy or sell 
one asset at price pif) which evolves as a function of time assumed to be 
discrete and measured in units of the time step At. In the simplest version of 
the model, each agent can either buy or sell only one unit of the asset. This 
is quantified by the buy state Si = +1 or the sell state Sj = —1. Each agent 
can trade at time t — 1 at the price pit — 1) based on all previous information 
including that at t — 1. We assume that the asset price variation is determined 
by the following equation 



a is the price volatility per unit time and r](t) is a white Gaussian noise with 
unit variance. The first term in the r.h.s. of (20) is the systematic price drift 
resulting from the possible imbalance between buyers and sellers. The impact 
function F(x) is such that F(0) = and is monotonically increasing with 
its argument as shown by recent empirical studies (Kempf and Korn, 1999; 
Plerou, et al., 2002; Lillo, et al., 2002; Bouchaud and Potters, 2003): perfect 
balance between buyers and sellers does not move the price; a larger (resp. 
smaller) number of buyers than sellers drive the price up (resp. down). An often 
used dependence is simply a linear relationship F(x) = fix (Beja and Goldman, 
1980; Bouchaud and Cont, 1998; Farmer, 2002). The second stochastic term of 
the r.h.s. of (20) accounts for noisy sources of price fluctuations. Taken alone, 
it would give the usual log-normal random walk process 

At time t — 1, just when the price pit — 1) has been announced, the trader i 




(20) 
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defines her strategy Si(t — 1) that she will hold from t — 1 to t, thus realizing 
the profit/loss (p(t)— p(t — l))si(t — 1). To define Si(t — 1), the trader calculates 
her expected profit Pg, given the past information and her position, and then 
chooses Si(t — 1) such that Pg is maximum. Within the rational expectation 
model, all traders have full knowledge of the fundamental equation (20) of 
their financial world. However, they cannot poll the positions Sj of all other 
traders which will determine the price drift according to (20). The next best 
thing that trader i can do is to poll her N(i) "neighbors" and construct her 
prediction for the price drift from this information. Note that, in this approach, 
the "neighbors" are by definition those who are polled by the trader according 
to her network of acquaintance. The trader needs an additional information, 
namely the a priori probability P + and P_ for each trader to buy or sell. The 
probabilities P + and P are the only information that she can use for all the 
traders who are not polled directly. From this, she can form her expectation 
of the price change. The simplest case corresponds to a neutral market where 
P + = P_ = 1/2. The trader % expects the following relative price change 



where the index j runs over the neighborhood of agent i. Notice that the sum is 
now restricted to the N(i) neighbors of trader %. The contribution of all other 
traders, whom she cannot poll directly, is one contribution to the stochastic 
term a rj(t). The restricted sum over the neighbors is represented by the star 
symbol. Her expected profit is thus 



The equation recovers (6) (without the term ai(t)G(t)) by identifying with 
Kij = K taken uniform. The evolution of opinions given by (23) is nothing 
but the dynamical version of the Ising model in which the sentiments {sj} are 
called "magnetic spins." Recall that the Ising model exhibits a phase transition 
between two phases: 

(1) For weak coupling between spins (or large idiosyncratic noise a), the spins 
take random signs and there is no majority opinion. The average opinion 
is zero. 

(2) Above a threshold K c for the average coupling coefficient (or below a 
threshold a c ), the spins align spontaneously along a preferred direction; 
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(22) 



The strategy that maximizes her profit is 




(23) 
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there is a non-zero majority opinion, which can be either +1 or —1 de- 
pending on the history. The properties of this transition and the sponta- 
neous symmetry breaking in the context of social agents is described in 
details by Sornette (2003, see Ch.5). 

In the extension (1) with (5), the coupling coefficients Kij between "spins" 
is allowed to change with time. During their dynamical evolution, one can 
expect that the coupling coefficients between agents' sentiments explore 
a large set of possible values according to a process similar to the increments 
of a random walk-like trajectory. At some times, a significant fraction of the 
may approach or pass above the critical value K c : this will lead, accord- 
ing to the mechanism of the phase transition in the Ising model, to the oc- 
currence of a majority of opinion and thus to strong herding behavior. The 
introduction of the dynamics of the coupling coefficient K^t) is reminiscent 
of the mechanism discussed by Stauffer and Sornette (1999) in the similar 
context of a percolation model of cluster of agents impacting the price evolu- 
tion Cont and Bouchaud (2000), in which the connectivity parameter is varied 
randomly in time. It is also providing a specific mechanism for the occurrence 
of critical times as explained by Sornette (2003, see Ch.5). The additional ex- 
istence of external news corresponds in the language of the Ising model to an 
external forcing "field" which help the opinion of the majority in the herding 
phase to bifurcate towards one or the other values ±1. However, the dynamics 
of (1) with (5) is clearly more complicated than for the Ising model due to the 
feedback of the external news and the collective decision process captured by 
the term (3r(t)G(t) on the coefficients of influences K^. 
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Fig. 1. Density distribution of returns T\ for a realization of the artificial stock 
market model formulated by Goncalves (2003) generated using 6 max = 0.22 ~ 0.24, 
Cmax = 0.14 ~ 0.15 and CV = 0.8 ~ 0.9 as recommended by this author. The time 
series of returns have been kindly provided by Gongalves. Our own simulations 
reproduce the same results. 
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Fig. 2. A typical example of the multimodal distribution for b n 
o- m ax = 0.045, and CV = 0.1. 
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Fig. 3. A realization of the logarithm of the price over 10 5 time steps generated 
using a = 0.2, 6 max = 0.3, er max = 0.03 and CV = 0.1 of the generalized artificial 
stock market model defined by (1), (4) and (10). 
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Fig. 4. Time series of the log-returns of the price shown in Fig. 3. 
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Fig. 5. (Color online) Empirical (solid lines) and theoretical (dashed thin lines) 
probability distribution density (in logarithmic scales) of log-returns at different 
time scales r of the price time series shown in Fig. 3. The log-returns r T are normal- 
ized by their corresponding standard deviations a T . The pdf curves are translated 
vertically for clarity. The thick dashed line is the Gaussian pdf. 
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Fig. 6. Autocorrelation function of the log-returns of the realization shown in Fig. 3. 
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Fig. 7. Autocorrelation function of the absolute value of log-returns of the realization 
shown in Fig. 3. The top panel show the correlation in linear-linear scale. The bottom 
panel plots the correlation function as a function of the logarithm of the time lag, 
as suggested by the multifractal random walk model (see text). 
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Fig. 8. The impact of a on the auto-correlation of the absolute values of the returns 
and of the returns. 
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Fig. 9. Scaling of the autocorrelation functions of |r r (t)| for different time scales r 
of the realization shown in Fig. 3. 
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Fig. 10. Scaling of the structure functions M q {r) of log-returns shown in Fig. 4 at 
different scales r for different orders q. 
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Fig. 11. Dependence of the scaling exponents £ q defined in (14) as a function of the 
order q of the structure functions M q (r) ~ r^ q . The concavity of £ 9 as a function q 
is the hallmark of multifractality. 
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Fig. 12. Average normalized conditional volatility (?\ t {t) / E[a 2 ] as a function of the 
time t — t s from the local burst of volatility at time t s for different log-amplitudes 
s in double logarithmic coordinates. 
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Fig. 13. Exponent ct(s) of the conditional volatility response as a function of the 
endogenous shock amplitude S for At = 1,2,4, and 8. 
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Fig. 14. Relaxation of superposed excess volatility after exogenous shocks obtained 
by imposing a very large news G(t s ) for At = 1. 
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Fig. 15. Five price trajectories showing bubbles preceding crashes that occur at the 
shifted time 0. The five time series have been translated so that the time of their 
crash is placed at the origin t = 0. 
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